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Stochastic master equation for a probed system in a cavity 
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We present a detailed derivation of the stochastic master equation determining the time evolution 
of the state of a general quantum system, which is placed inside a cavity and subjected to indirect 
measurements by monitoring the state of electromagnetic radiation transmitted through the cavity. 
The derivation is based on the physics involved and the final result is stated in terms of the phys- 
ical parameters of the setup. To illustrate the predictions contained in the equation, we solve it 
analytically for a specific system, and we demonstrate quantum jumps and freezing of the internal 
coherent dynamics of the system as a result of continued measurements. 
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I. INTRODUCTION 

The time evolution of the state of an isolated quan- 
tum mechanical system is governed by the Schrodinger 
equation. If the system is subjected to a measurement, 
the state of the system collapses onto an eigenstate of 
the measured observable, and the measurement outcome 
is the corresponding eigenvalue. Instead of performing a 
measurement directly on the system, one may also allow 
the system to interact with an auxiliary system, which 
is then subjected to a measurement. This measurement 
procedure reveals only partial information on the state 
of the original system and is very useful in protocols to 
prepare a system in a certain quantum mechanical state 
and to achieve quantum nondemolition measurements of 
an observable. The system could for instance be clouds 
of atoms or few atoms that are probed with a beam of 
light [H, 0) [1, 0| or an electromagnetic field mode in a 
cavity probed with a beam of atoms [H . 

When a system interacts with the surroundings, it is 
not practical to keep track of the complete state of both 
system and surroundings, and instead the degrees of free- 
dom of the surroundings are traced out. In this case the 
state of the system is no longer pure and must be de- 
scribed by a density operator, whose time evolution is 
determined by the master equation. For a system sub- 
jected to indirect measurements, the time evolution de- 
pends on the actual outcomes of the measurements, and 
since these are probabilistic in nature, they are handled 
mathematically by introducing stochastic variables into 
the master equation. This allows one to calculate both 
the state of the system conditioned on a given set of mea- 
surement outcomes (a given realization of the stochastic 
variables) and the probability to observe that particular 
sequence of measurement results @. 

The stochastic master equation is often derived in very 
general settings, and it has been shown that all master 
equations must be on the so-called Lindblad form 0] in 
order to preserve complete positivity of the density op- 
erator. The generality is, however, achieved at the cost 
of introduction of abstract measurement operators and 
measurement strengths. In contrast, in the present pa- 



per we assume a concrete measurement scheme and de- 
rive the stochastic master equation directly from known 
physical interactions. The system under consideration is 
placed inside a cavity and probed with electromagnetic 
radiation by sending photons into the cavity from one 
side and performing homodyne measurements on the field 
leaving the cavity on the opposite side. The purpose of 
the cavity is to reflect the probe light several times before 
it is detected whereby the effective interaction strength 
between the light field and the system is increased as we 
shall demonstrate below. 

Probing of atomic systems by their effect on the trans- 
mission properties of an optical cavity have a long his- 
tory in quantum optics. Early studies focused on the 
field-atom dynamics, leading, for instance, to the normal- 
mode splitting of the transmission resonance Q and on 
photon statistics 0], while more recent work has shown 
the possibility to observe the spatial motion of individ- 
ual atoms trapped inside the cavity by the probing beam 
itself [l(| [HI • See also the work on optically transported 
or guided atoms [H, [H, [3] . More recently, optical cavi- 
ties have been introduced in experiments to probe Bose- 
Einstein condensates [HI, [l7j ■ 

A stochastic master equation for a setup involving a 
cavity and a homodyne detector has been derived using 
a rather different approach in Ref. (l8| . In that paper, 
however, the aim was not to use an auxiliary system to 
perform indirect measurements but to determine the time 
evolution of the field in a cavity, when the light leaking 
out of the cavity is subjected to measurements, and thus 
the probing light and the probed system, which are cru- 
cial ingredients in the present paper, were not included 
in the analysis. 

The paper is structured as follows. In Sec. |TT] we de- 
scribe the physical setup under consideration and discuss 
the time evolution of the density matrix. If the dynam- 
ics is slow compared to the round trip time of light in 
the cavity, we can use a continuous description and de- 
rive the stochastic master equation in Sec. IIIII In Sec. 
IIVI we comment further on the derived equation, and in 
Sec.[V]we apply it to specific systems to illustrate explic- 
itly how the state of the system is gradually collapsed 
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by the measurements and to investigate how the relative 
strength of the probing and the coherent dynamics of the 
system influences the time evolution. Sec. I VII concludes 
the paper. 



II. MODEL OF THE PROBING PROCEDURE 

The probing procedure applied is shown schematically 
in Fig. [TJ The system is placed inside a four sided ring 
cavity. The probe light enters the cavity from the left, it 
interacts with the system, and eventually it leaks out of 
the cavity, where the transmitted light is subjected to a 
balanced homodyne measurement. In order to describe 
the quantum state of the light field we divide all the light 
beams into small segments of duration dt as illustrated 
in Fig. [2] and treat each segment as a single mode. This 
is valid provided dt is chosen sufficiently small. The state 
of the light and the system is then specified as the col- 
lective state of all the light modes and the system. To 
avoid keeping track of an enormous number of modes, we 
assume that the probe light is in a coherent state before 
it enters into the cavity. The modes of the probe beam 
that have not yet reached the cavity are then in a prod- 
uct state, and their time evolution is independent of the 
dynamics inside the cavity and the outcome of the ho- 
modyne measurements. These modes can consequently 
be excluded from the analysis until the time, when they 
arrive at the cavity. The vacuum field incident on the 
beam splitter BS2 in Fig. [JJ and the coherent state local 
oscillator field are also product states, and again it is suf- 
ficient to consider the modes that are currently entering. 
The modes that are detected in the homodyne detector 
are traced out after the measurements, and the unob- 
served modes that leave the cavity at BSi are also traced 
out in each time step. The number of active modes is 
then constant in time, and in the following we denote 
the density operator of the light modes inside the cavity, 
the light modes between the cavity and the detectors, 
and the system at time t by p(t). We note that the anal- 
ysis can also be carried out for a squeezed input field by 
including an optical parametric oscillator in front of the 
cavity in Fig. [TJ For a suitable theoretical description of 
the optical parametric oscillator see Ref. [l9[. 

To determine the time evolution of p(t) we express 
p(t + dt) in terms of p(t). The relevant modes are the 
modes included in p(t) and the mode of the probe field, 
the vacuum state mode, and the local oscillator mode 
that reach beam splitters BSi, BS2, and BS3, respec- 
tively, in the time interval between t and t + dt, and the 
state of these modes is p(t)&) pi n (t)®\0) (0\®\a){a\. In the 
time between t and t + dt several interactions take place. 
The state of the two modes that hit the beam splitter 
BSi are mixed, and this process is described by the uni- 
tary operator U\. Similar transformations, described by 
the operators U2 and U3, occur at beam splitters BS2 and 
BS3. The system itself undergoes evolution during dt and 
the light mode that passes the system at time t interacts 
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FIG. 1: Model of the setup. The cavity consists of the two 
mirrors Mi and M2 and the two beam splitters BSi and BS2. 
The probe light pi n (t) enters the cavity through BSi. The 
cavity field interacts with the system and leaks out of the 
cavity through BSi and BS2 (|0) is a vacuum state). The 
light transmitted through the cavity is subjected to balanced 
homodyne detection. BS3 is a 50:50 beam splitter, \a) is a 
strong local oscillator, and the two detectors register n and m 
photons, respectively. The measurement readout is k — n—m. 



dt 

i > 



FIG. 2: The light beams are divided into temporal modes of 
(infinitesimal) duration dt. 



with it. The corresponding infinitesimal time evolution 
operator is Uh = 1 — iHdt/h, where H = H sys + Hi is 
the sum of the system Hamiltonian H sys and the Hamil- 
tonian Hi for the interaction between the system and the 
light field. We note that the system could be subjected 
to manipulations that depend on the state of the system 
and light field and such feedback terms would be included 
in H sys . In the following we disregard decay of the sys- 
tem, since it does not add any interesting points to the 
analysis, and it is easily included by adding a term to the 
final equation (see for instance [|| for a derivation of the 
relevant term for a system consisting of a single two-level 
atom). If the system is an extended object that inter- 
acts strongly and, for example, depletes the light field, 
it would be relevant to slice the system into small pieces 
and consider interactions between the system and several 
light modes in each time step, but we ignore such compli- 
cations here. We shall also assume that the mirrors Mi 
and M2 induce phase shifts of tt/2, which are taken into 
account through the operators Um-l and Um 2 - Finally, 
the two modes that hit the detectors in the time inter- 
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val dt are projected on photon number states \n) and 
|m), respectively. Even though the detectors are able 
to resolve the exact number of photons microscopically, 
the only macroscopically available measurement readout 
is the difference k = n — m between the number of de- 
tected photons in the two detectors. We thus have to 
average over all possible values of n and m that lead to 
the observed value of k. Putting all the transformations 
together, we obtain the density operator at time t + dt 
conditioned on the measurement of a difference of k pho- 
tons in the interval from t to t + dt 

p(t + dt) = ■p-Tn(22{m + k\{m\UsUiU Ma 

m 

U Ml U 2 U H p(t)®Pin(t) ® |0)(0| <8> \a)(a\ 

UlUp^ Ul 2 UMW \m + k)), (1) 

where Pk is the probability, determined from the normal- 
ization of p{t + dt) , to obtain the measurement outcome 
k, and TVi denotes the trace over the unobserved mode 
leaving the cavity at beam splitter BSi. 

The time evolution of p(t) is completely specified by 
Eq. ([T]) , and in principle one can start with a given initial 
state and iterate {1} on a computer. This task does, 
however, soon become unwieldy due to the large number 
of modes that are involved. In the next section we show 
how Eq. JT]) can be simplified considerably and rewritten 
as a stochastic differential equation by invoking a few 
assumptions. 



III. DERIVATION OF THE STOCHASTIC 
MASTER EQUATION 

The key assumption in the following derivation is that 
all changes of the state of the light field takes place on a 
time scale that is large compared to the round trip time 
t of light in the cavity, which is valid if the cavity is suffi- 
ciently small. In this case the temporal width of the light 
modes can be chosen as large as r, and, in particular, the 
field inside the cavity can be treated as a single mode. 
The assumption requires that the coupling between the 
system and the light field is not too large, since, for in- 
stance, the fraction of the light absorbed or emitted dur- 
ing a single round trip must be small. The finesse of the 
cavity must also be high and the total number of pho- 
tons in a segment of duration r of the input beam must 
be small in order to avoid that a significant number of 
photons leak out of or enter into the cavity within a time 
t. Finally, variations in, for instance, the input field or 
the interactions between the system and the light field 
that take place on a time scale small compared to r are 
not allowed, since they will be smeared out. Actually, the 
input field typically has such fast variations because the 
size of the cavity is large compared to the wavelength of 
the input light, but, provided the frequency of the radia- 
tion is not too far from a cavity resonance, this problem 



is easily circumvented by moving into a frame rotating 
with the relevant resonance frequency of the cavity, and 
in the following all frequencies are to be measured rel- 
ative to the cavity resonance. With these requirements 
we can treat r as an infinitesimal quantity, and below we 
shall denote r by dt. For simplicity, we shall also assume 
that the homodyne detector is placed within a distance 
ct from BS2, where c is the speed of light, since in that 
case pit) includes only the state of the system and a single 
mode of the light field (the cavity mode) . 

When r is infinitesimal, the transmission t\ (t|) of 
beam splitter BSi (BS2) must also be infinitesimal, and 
we define k, by 

t1 = Kidt, i = 1,2. (2) 

This allows us to rewrite the unitary operator 

TT jlLfAa jJLh^h — i tan -1 ( — ) f a^b+ab^) /oX 

Ui = e 2 e 2 e W/V ^ (3) 

representing BSi, as 

^1 - (atS + aft) - h\id r b + dftf^j (4) 

to first order in dt, where r\ is the reflectivity of BSi, a 
is the field annihilation operator of the cavity field, and 
b is the field annihilation operator of the mode of width 
dt of the input field that arrives at the beam splitter at 
time t. Note that has been chosen such that after the 
interaction, the a-mode is still the cavity mode, while the 
6-mode is the mode leaving the cavity. Similarly, we have 

^1 - it 2 ia)c + ac f ) - ^tUd^c + ac f ) 2 ^j , (5) 

where c is the field annihilation operator of the mode 
leaving the cavity at BS2. The operators Um^ and Um 2 
both act on the cavity mode and thus 

U Ml = U Al2 = e^ ata . (6) 

Since the number of photons in a segment of duration r 
of the input beam is assumed to be much smaller than 
one, we can write 

ftn(t) - c in , 00 |0)(0| +C in ,i O |l)(0| + C in ,oi|0)(l| 

+ Cin,20 |2)(0| + Cin,02 |0)(2| (7) 

to first order in dt, where C; n) oo = 1 — Cm.n is of or- 
der unity, Cj rii io = c* n 01 is of order \fd~i, and Ci n ,n and 
Cin,20 = c* n 02 are of order dt. 
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Putting the relation Uh = 1 — iHdt/h and Eqs. ([4]), 
©, ©, and into Eq. (JIJ, we find 

1 / i 
p(i + dt) = — u 0Q p(t) - u 00 -[H,p(£)]<ii + uoo^i 
n V " 

([at,p(i)]Tr(6p in (t)) - [a,p(t)]Tr( Pm (t)6t)) 
+ u w t 2 ap(t) + uoit2p(t)a* 
+ ^"oo^i {—a)ap(t) — p{t)a) a + 2ap(t)a^) 

+ ( — ttoo^<zp(t) — itoop(£)a^& + 2unap(t)a^ 

+ \/2u 20 a 2 p(t) + V2u 02 p(t)(a t ) 2 )^ , (8) 

where 

u pq ee ^(m + fc|(TO|[/3|a)|p)< 9 |<a|C/ 3 t |m>k + ^) (9) 

m 

and |p) and |g) are photon number states. For a strong 
local oscillator (|a| 2 3> 1) we can apply the approxima- 
tion 



1 n -a 1 '-^ 2 
n! \/2ttp 



(10) 



and turn the sum in Eq. ([5]) into an integral, which leads 
to (a = He 1 *) 



"10 



k 2 \ 



-ifce-^ 1 / fc 2 \ 

V2p 6XP V 4p~/ 



U20 = - 



(k 2 - 2/i)e" 2 ^ 1 

f= „ ■ exp , 

2a/2^ 20F7I V V 



k 2 1 / fc 2 ^ 
271 2^/^71 6XP V <W ' 



(11) 
(12) 

(13) 
(14) 



From the normalization of Eq. (|8|) we find the probability 
to obtain the measurement outcome k 

1 f k 2 \ f ike-** , , „ 



ifce 



1 



i 2 Tr (p(t)at) + -f 



2// 



1 ) f2Tr(o/3(i)t 



-2i« 



: ^Tr (a 2 p(t)) - e 2l0 Tr (p^a 1 ) 2 ) 
Comparing this to 

(k-S) 2 



(15) 



-. exp 



4(M + e) 



20?M eXP ( 4p) ( 1+ 2 A1 + (2 Ai 2/J 



for <5 <C fc and e <C p, it is apparent that to first order 
in dt, is a Gaussian distribution. S is of order y/di, and e 
is of order dt, so to order s/dit we can replace fc/ ' \/2p by 

-A= = -ze-^ 2 Tr (ap(t)) + ze ,; ^ 2 Tr (p(i)a f ) + 
2p Vc?i 



(17) 

where dW is a stochastic variable that has a Gaussian 
probability density distribution with zero mean value and 
variance dt, and we have put a hat on k to emphasize 
that k is now to be regarded as a stochastic variable that 
assumes the value k with probability From the Ito 
calculus rule dW 2 = dt (see 0) it furthermore follows 
that 



k 2 = 2(i 



(18) 



to zeroth order in dt. Inserting Eqs. ©, (fTTj). fl2|. (fT3|) . 
d, (USD, dUD and 



Tr 



(Sp in (t)) = Tr (pinW^) * = c ln , 10 ee pVdt (19) 



into Eq. (|8|) we finally obtain the stochastic master equa- 
tion in the form 



p(t + dt)=p(t)--[H,p{t)]dt 

+ ^fk~[[ci\ p(t)](3dt - V^TIo, p{t)](3*dt 
- ie'^y/^ (ap(t) - Tr (ap{t)) p{i)) dW 
+ ie*fy*a (p(*)a f - Tr (p(t)tf) p(t)) dW 
+ - (k x + K 2 ) (-a f ap(t) - p(i)a f a + 2ap(t)a t ) dt. (20) 

The second term on the right hand side represents the 
evolution due to the system Hamiltonian and the interac- 
tion between the system and the radiation, the third and 
fourth terms arise due to the feeding of probe light into 
the cavity at BSx, the fifth and sixth terms include the 
knowledge obtained from the homodyne measurements, 
and the seventh term describes the decay of the cavity 
mode due to transmission through BSx and BS 2 . We 
note that the approximations leading to Eq. ([20]) can be 
stated more precisely as \a\ 2 3> 1, kit <C 1, k 2 t <C 1, 
y/Ki\/3\r < 1, and \(ipi\H\<ip 2 )\T/h < 1, where \tpi) and 
\ip 2 ) represent arbitrary state kets. It is allowed that j3 
and H arc time dependent, but the variation within a 
time interval of length r should be small. 



IV. FURTHER REMARKS 

So far we have assumed a lossless setup, but it is easy 
to incorporate effects of losses. To account for losses in 
the light field, we only need to replace the perfect mirror 
M2 by a partially transmitting beam splitter and include 
a beam splitter in front of the ideal homodyne detector 
with transmissivity 77 = r\Df]p, where rju is the efficiency 
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of the detector and 1 — rjp is the propagation loss between 
the cavity and the detector, and in this case Eq. ([ST)]) 
generalizes to 



p(t + dt) = p(t)--[H,p(t)]dt 

+ ^[a\p(t)] 0dt-y/Ki[a,p(t)] (3*dt 



- ie-^y/rpvi (ap(t) - Tr (ap(t)) p{t)) dW 
+ ie^y/rjRi (p(t)a^ - Tr (p(t)a f ) p{t)) dW 

+ —- k (-a)ap(t) - p{t)a)d + 2ap(t)a)) dt, (21) 

where k = Ki + n 2 + Kl, and kl, defined in analogy to 
Eq. @, is the decay rate due to intra cavity loss. It is 
also possible to account for loss and decoherence of the 
quantum system in the cavity by introducing appropriate 
damping terms in Eq. (I2ip . 

If there is no decay of the system and no cavity loss 
and all light emerging from inside the cavity is detected 
with unit efficiency detectors, the dynamics will preserve 
the purity of an initially pure state of the atoms and the 
cavity field, and the stochastic master equation can be 
rewritten as a stochastic Schrodinger equation. This is 
convenient since it is significantly easier to propagate a 
wave function in time than a density operator. The situ- 
ation can be achieved with the setup in Fig. [1] by taking 
the limit K\ — > and |/3| — > oo while ^/k7|/3|t is kept 
small compared to unity. Alternatively one could mea- 
sure both the light reflected and transmitted from the 
cavity or replace BS2 by a perfectly reflecting mirror and 
subject the light leaving the cavity at BSi to homodyne 
detection. In the latter case the stochastic master equa- 
tion is given by Eq. (I20[) with k 2 and K\ + Ki replaced 
by Ki and cf> replaced by <fi + and the stochastic 
Schrodinger equation reads 

\ip(t + dt)) = (l - ^Hdt + fiy/Kia^dt - f3*y/K~[ddt 
- e'^y/^ia - (a))dW 
- y (a* a - 2a{a)) + {a){a*))dt} |V>(t)>, (22) 

where (•) denotes expectation value. 

Equation (|21[) is a nonlinear equation since Tr (ap(i)) 
and Tr (p(t)a>) depend on the state pit). In general, H 
may also depend on pit) (for instance if a state dependent 
feedback is applied), but in the special case of a state 
independent Hamiltonian, it is possible to transform Eq. 
(I2ip into a linear equation by application of the method 
presented in Ref. [y]. The nonlinear terms in Eq. (|2"Tj) 
appear due to the first two terms on the right hand side 
of Eq. (fT?]) and due to the normalization of the state. If 
we simply remove the factor 1/Pfc in Eq. (8]), the norm 
of pit + dt) is the probability to obtain the measurement 
outcome k for the time interval between t and t + dt. 
Furthermore, the right hand side of the master equation 
is determined completely by the value k assumed by k 



(or equivalently by the value dW assumed by dW) and is 
thus not changed if the probability density distribution of 
dW is changed. Only the probability to obtain the value 
of dW that leads to a specific value of k is changed. We 
can thus change the probability density distribution of 
dW provided we accept that the probability to obtain 
the state pit + dt) at time t + dt is the norm of the state 
and not the probability to obtain the required value of 
dW . We exploit this freedom to define the new stochastic 
variable dy by 



k 



dy_ 

dt 



(23) 



and assume that the probability density distribution of 
dy is a Gaussian distribution with zero mean value and 
variance dt. If we omit the factor 1/Pfe in Eq. ([5]) and 
insert Eqs. (JTTJ, (12]), (JTSJ), (O, and (23]), we obtain a 
linear equation in pit). The trace of pit + dt) integrated 
over all possible realizations of k is unity, but since k 
is expressed in terms of dy, we would like the trace of 
p(t + dt) integrated over all possible realizations of dy 
to be unity, and we thus multiply the right hand side by 
\j2pjdt. Finally, to ensure that p{t + dt) approaches pit) 
in the limit, where dt approaches zero, we divide the right 
hand side by exp (— fc 2 /(4/i)) /y / 2ndt, which is precisely 
the Gaussian probability density Pd y for dy. With these 
changes Eq. (2~Tj) reduces to 



pit + dt) = pit)--[H,p[t)]dt 

+ v^I [a*, pit)) (3dt - ^ [a, pit)} 13* dt 

- iy/v^ (ap(i)e -40 - p(t)a f e # ) dy 
+ (-a f ap(t) - pit)a)d + 2ap(t)a t ) dt. (24) 

Due to the division by Pd y , the probability density to 
obtain a given state p{t + dt) at time t + dt is now the 
product of the norm of pit + dt) and the probability den- 
sity Pdy to obtain the required value of dy. Since the 
stochastic variables corresponding to different time steps 
are independent and since Eq. (|24p is linear in pit), the 
probability density to obtain the specific state after N 
time steps is simply the probability fL—i ^dvt t° obtain 
the realizations dyi , dy 2 , ■ ■ ■ , dyN of the stochastic vari- 
ables multiplied by the norm of the state obtained from 
this realization summed over all realizations that lead to 
the desired state. In particular, if the state after N time 
steps only depends on the sum y = ^ i=1 dyi, then the 
probability density to obtain this state is 



P = 



1 



y/2irNdt 



exp 



V 



2Ndt 



Tr(p(t + Ndt)). (25) 



In the preceding section we have assumed that the 
transmitted light is observed by a homodyne detec- 
tor, but a similar derivation may be carried out for an 
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avalanche photo diode detector, which projects the in- 
finitesimal modes of the transmitted light on either vac- 
uum or a single photon state, and for completeness we 
state the result 



p(t + dt) =p(t) - -[H,p(t)]dt 



ni \a* , p(t)] (3dt — \f~k~i [a, p(t)] f3*dt 
— rjK2{ap{t)a} — Tr (a p(t) a* )p(t))dt 

+ -n(-a f ap(t) - p{t)a)a + 2ap(t)a f )dt. (26) 

dN is a stochastic variable, which assumes the value 1 
with probability r]K2TT(ap(t)a' ! )dt and the value with 
probability 1 — r]K2H(ap(t)a^)dt. 



V. APPLICATIONS 
A. The empty cavity 

We first consider the time evolution of the state of 
the light field in the absence of any quantum system 
in the cavity. This situation leads to a coherent state 
cavity field because beam splitters transform coherent 
states into coherent states and because homodyne mea- 
surements also preserve the coherent state nature of the 
cavity field (see e.g. Ref. [20] ) - Inserting H = and 
p{t) = C(t)\Z(t))(£(t)\ in Eq. flST)), we obtain 



dt ~ 2 U ] 



and if (3 is time independent, which implies that the input 
light is on resonance with the cavity, £ (t) approaches £ = 
2- N /Ki/3//t for t > 2/k. 

We note that £ is a factor of 2^fk~\j [yfrn) larger than 
the coherent state amplitude y/r(3 of a segment of the in- 
put field of length ct, and, in the limit of very weak cou- 
pling between the light field and a system, the presence 
of the cavity enhances the Rabi frequency of transitions 
between different states of the system by the same fac- 
tor. For K3 = and k% = K2 the factor reduces to l/t2, 
which equals the square root of twice the average num- 
ber of round trips of a photon in the cavity in absence of 
the system and in absence of the input field. The factor 
of two appears due to destructive interference between 
the input field and the cavity field for £(t) = 2y/K\/3/K, 
which ensures that no light is lost at beam splitter BSi 
in the presence of the input light. 



B. Analytical solution for a simple system 

As a nontrivial application of the stochastic master 
equation we analyze the setup proposed in Ref. [H to gen- 



erate Dicke states. Here the system consists of N nonin- 
teracting identical atoms each with two ground state lev- 
els I/) and \g). The cavity field couples the level |/) to an 
exited level |e), and below it is assumed that the radia- 
tion is sufficiently off-resonant to avoid population of the 
exited level. In that case the interaction merely shifts the 
phase of the light field by an amount that is proportional 
to the number of atoms in the state |/) , and since the ho- 
modyne detector is sensitive to the phase shift, the state 
of the monitored system slowly approaches an eigenstate 
of the operator h, which counts the number of atoms in 
the state |/). For atoms all coupling with equal strengths 
to the light field, the relevant Hamiltonian reads 



H = hgaJah, 



(28) 



where g represents the strength of the coupling. We as- 
sume that the atoms are initially prepared in a state, 
which is symmetric under exchange of any two atoms, 
and we denote the symmetric state with n atoms in the 
state I/) by \n). Since the Hamiltonian (|28p preserves the 
symmetry under exchange of any two atoms, this allows 
us to use the restricted basis consisting of the states |n) 
with n = 0, 1, N. 

To simplify the stochastic master equation, we write 
the density operator of the system and the cavity mode 



as 



N N 



(29) 



where \n) and |m) refer to the state of the system and 
Pnm are time dependent cavity mode operators. Inserting 
this into Eq. ([24]) , we obtain the (TV + l) 2 independent 

linpfir prmfltinnQ 



(27) linear equations 



dp r , 



dt 



it; 



= -ig(na> ' dp nm - mp nm a [ <i 

-i(f> 



iy/r]K 2 {ap nm e 
1 



Ki [a,p nm \ (3* 
d£ 
dt 



o aV' 



-K I 



p nm a)d + 2ap nm a t ) . (30) 



To determine (n|p sys |n), where p sys is the normalized 
density operator for the system, it is sufficient to solve 
Eq. (f3"0"l) for n — m. If the initial state of the system is 
|n)(n|, we can replace the operator h in the Hamiltonian 
(|28| by the number n, whereby the system is effectively 
reduced to a phase shifter, and it follows that the cavity 
field is in a coherent state. Equation (|30[) for n = m 
is, however, not mathematically different in the general 
case, and it is thus solved by 



Pnn=C n (t)\Ut))(Ut)l 

Inserting Eq. (JHTJ) into Eq. $5U§, we find 

d£n(*) 



(31) 



dt 



(| + ign) £ n (t) + (32) 
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and 



dC n (t) = ^jJT 2 (-te-^„(t) + ie<*£{t)) C n (t)dy, (33) 
and thus, for a time independent /3, 



u*) = U»(°)- 



k + 2i<?n 



and 



'(34) 



C„(t) = C„(0)exp 



2 7o 



r n (ir*' , (35) 



where r n {i) 



Hn{t) 



ie 



i(f> 



£(*))■ 



The general solution simplifies considerably if the de- 
tector is turned off (77 = 0) until £, n (t) has obtained its 
steady state value 



2y r Ki/3 1 — 2ign/ k 
k 1 + Ag 2 n 2 /n 2 



(36) 



since in that case r n (t') = r n can be moved outside the 
integrals, and C n (t) depends only on the sum y of all the 
infinitesimal increments dy. (n\p sys \n) does not change 
as long as 77 = 0, and below we thus simply assume that 
£„(0) is the steady state value. We note that this shift 
is insignificant if the changes of the state of the system 
induced by the homodyne measurements take place on 
a time scale that is much longer than which is of- 
ten the case, because the coupling strength g of the off- 
resonant interaction between the system and the cavity 
mode is typically very small. 




-200 -150 -100 -50 

ry 



FIG. 3: (a)-(c): Probability density P to measure the inte- 
grated signal y in the interval from to t as a function of y 
for r 2 t = 1, 10, and 50, respectively, (d)-(f): {n\p sys \n) for 
n = 0, 1, 2, 3, 4 as a function of y for the same values of r 2 t. 
The values of d(0) are those obtained for the initial state 
IV'sys) = ((I/) + ls))/v / 2)® 4 of the system, and it is assumed 
that 2g <C k. 



C. Probing as a state preparation tool 

We next consider in more detail how the measurement 
of the phase shift of the transmitted light can be used to 
prepare different types of quantum mechanical states. 

1. Dicke states 

To prepare a Dicke state we choose = (correspond- 
ing to a measurement of the p-quadrature of the field), 
since it follows from Eq. (|36[) and the definition of r n 
that those measurements are most sensitive to the in- 
duced phase changes in the limit 2g/n <C I and j3 — (3*. 
In this case 

C„(t) = C„(0) exp(-r„y - r 2 n t/2), (37) 

with 



and, from Eq. (|25|) . the state for which y assumes the 
value y is obtained with probability density 

P = V%Eexpp + ^ 2 y (39) 

P and the state preparation fidelity (n\p sys \n) = 

CnW/Em=o^™W are plotted as functions of y at dif- 
ferent times for C 4 (0) = N\/i\/(N - i)\/2 N and N = 4 
in Fig. [3l neglecting the term 4g 2 n 2 in the denominator 
of Eq. (|38p . such that r n = rn, where r is independent 
of n. Within this approximation P consists of a sum 
of five Gaussians separated by rt and of width y/t, and 
the transition from overlapping Gaussians at small r 2 t to 
well separated Gaussians at large r 2 t is apparent in the 
figure. The plots of (n\p sys \n) illustrate how the state 
of the system is gradually collapsed onto an eigenstate 
of the operator n, and it is clear that each peak in P 
corresponds to a specific value of n if r 2 t is large. 
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2. Quantum super-position states 

It is also possible to use the setup to generate a su- 
perposition of two Dicke states as, for instance, the max- 
imally entangled state (|0) + e l6 \N))/y/2. To do so we 
include an additional phase shift inside the cavity such 
that the Hamiltonian reads 

H = hga)a(n — N/2). (40) 

With this Hamiltonian the states \n) and \N — n) are 
indistinguishable if the re-quadrature of the field is mea- 
sured (0 = — 7r/2 in Eq. J20|)), an d thus, for the initial 
state |VVs) = ((|/) + \g))/\/2)® N , a superposition of |0) 
and \N) is generated after sufficiently long measurement 
time with probability 1/2 N ~ 1 . To achieve a pure atomic 
state, it is required to work under lossless conditions, i.e., 
under conditions where Eq. ([22]) is valid, and it is also 
important to avoid entanglement between the cavity field 
and the atoms at the final time, since the cavity field is 
to be traced out to obtain the atomic state. The last re- 
quirement is easily fulfilled by turning off the probe field 
at some point, letting the cavity field decay to the vac- 
uum state while keeping the detector turned on. We note 
that the protocol is more demanding than the protocol 
to generate Dicke states due to the requirement of loss- 
less conditions and because the difference between r n and 
r n+ \ is second order in 2g/n 1 which leads to an increase 
in the required measurement time. 



D. Quantum Zeno effect and quantum jumps 

In the preceding subsections the state of the atoms was 
simply detected by an indirect continuous measurement, 
but we now add the complication of internal dynamics by 
allowing transitions between adjacent |n) states. Specif- 
ically, we consider the Hamiltonian 

N 

H = tigb^an + Hg s YJ(ct,- i+ + cr,,-), (41) 

i=l 

where cr^+ = \f)u(g\ and <t^_ = \g)u{f\, and i refers to 
atom number i, which, depending on the level structure 
of the atoms, may be realized by applying a microwave 
field propagating in a direction perpendicular to the cav- 
ity axis. Since 

N 

^2 *<,+ !«} = V(n + l)(N - n)\n + 1), (42) 

and 

N 

= Vn(N + l-n)\n-l), (43) 

i=l 
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FIG. 4: (0|p 3ys |0} as a function of scaled time g s t for N = 1, 
ki = K2 = 0.5 k, kl = 0, T) = 1, <j> = 0, g = 0.2 k, 
P = 0.2 v 7 ^, N p — 3, and initial state p(0) = |0}|i>ac)(mc|(0|, 
where \vac) is the vacuum state of the cavity mode. For weak 
driving with g 3 — 0.001 n (red curve) the continuous mea- 
surements inhibit the coherent evolution between the states 
|0) and jl), while for strong driving with g s = 0.05 k (green 
curve) the time evolution of the state almost follows the de- 
terministic evolution obtained for g — (dashed curve). The 
blue curve (g s = 0.005 k) illustrates the behavior of the state 
of the system in the intermediate regime. The same noise 
realization is used to compute all curves. 



we obtain the stochastic master equation for the system 



and the cavity field by adding the terms 

-ig s ^n{N + 1 - n)p n _i jm (44) 
-ig sy /(n+l)(N-n)p 

n+l,m 

(45) 

+ig s ^/m(N + 1 - m)p„, m _i (46) 
+ig s ^(m+ 1){N -m)p n , m+1 (47) 



to the right hand side of the master equation for the 
Hamiltonian in Eq. (f28|) . 

The added system Hamiltonian drives a coherent evo- 
lution between the different \n) states of the system, but 
if the measurements are sufficiently strong, the continu- 
ous back action of the measurements on the state of the 
system blocks the coherent evolution. To demonstrate 
this effect, which is known as the Quantum Zeno effect, 
we use the Milstein scheme (see Ref. [2l|) to solve the 
nonlinear stochastic master equation numerically for a 
particular realization of the measurement outcomes. The 
accuracy of the integration may be checked for g s = 
by application of the analytical results obtained in the 
previous subsection. We use the Fock state basis for 
the cavity mode and neglect terms in the density op- 
erator with more than N p photons in the cavity, where 
N p > 4 Ki \P\ 2 /k 2 (see Sec. lVAj) . The results are shown 
in Fig. SI Since (l|p sys |l) is not exactly zero, there is a 
small probability that the system switches to state |1), 
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FIG. 5: Quantum jumps for a single atom in a cavity. The 
parameters are the same as in Fig. \%\(g s ~ 0.001 «). 

where it is subsequently stabilized, and if the dynamics 
is integrated for sufficiently long time, several such quan- 
tum jumps between |0) and |1) will occur (Fig. [5]). 

We note that for general N the states with n = or 
n = N are easier to stabilize by measurements, since the 
square root factors in Eqs. (gH), (j45]), (j46]), and ((47j) are 
smaller for n close to or N than they are for n close to 
N/2. 

Finally, as a simple example of feedback, we note that 
it is possible to generate a particular eigenstate of n by 
increasing g s whenever (n\p sys \n) becomes small and de- 
creasing g s to zero whenever (n\p sys \n) becomes large. 



VI. CONCLUSION 

In conclusion we have considered the influence of a spe- 
cific indirect continuous measurement on the state of an 
abstract system, and we have derived the master equa- 
tion determining the time evolution of the state of the 
system conditioned on the measurement outcomes. This 
equation provides a tool to analyze a wide variety of sys- 
tems and phenomena in detail, and we have applied it to 
demonstrate the collapse of the state of a system onto an 
eigenstate of the measurement operator and to demon- 
strate quantum jumps and blocking of the internal dy- 
namics of a system due to continuous measurements. 

Indirect measurements are important because they 
constitute a special class of state manipulations that are 
very useful to prepare systems in specific quantum me- 
chanical states. The crucial point is that the complete 
collapse of the state due to back action of a measurement 
performed directly on the system is avoided by perform- 
ing instead the measurement on an auxiliary system that 
has interacted with the system. We showed examples 
with simple atomic level structures, but we emphasize 
that our analysis applies to quite general quantum sys- 
tems and, for example, to multi-level atoms with shelv- 
ing states leading to macroscopic quantum jumps in the 
transmitted intensity [12, [23j]. A further step towards 
achieving control on the quantum state of a system is 
to apply feedback [24], HE IHl by subjecting the system 
to disturbances that depend on the outcome of measure- 
ments. This situation may also be handled by the formal- 
ism developed in the present paper. The derived stochas- 
tic master equation is thus a valuable tool to investigate 
further state preparation protocols. 
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